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Abstract 

We explore the origin of the flat spectrum seen in some T Tauri stars by considering a three-component 
structure: a central star, a circumstellar disk, and a dusty halo. The radiative energy transport is faithfully 
treated by solving the angle- and frequency-dependent radiative transfer equation in two space dimensions 
assuming axisymmetry, and hence the radiative equilibrium temperature in the disk and halo is determined 
simultaneously. The disk is effectively heated by the scattering and reprocessing of stellar radiation through 
the halo. The large mid- to far- infrared excess originates from the photosphere of the warmed disk, resulting 
in a flat spectrum, as observed. The halo which we consider is observed as a compact reflection nebula, 
and is discriminated from extended, disk-like envelopes around flat-spectrum T Tauri stars. We show that 
the overall spectral shape of flat-spectrum T Tauri stars can be reproduced by the present disk-halo model. 
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1. Introduction 

It seems reasonable to suppose that the infrared ex- 
cesses of T Tauri stars can be attributed to thermal dust 
emission from circumstellar disks. By modeling the ob- 
served spectral energy distributions (SEDs), one can de- 
rive the disk properties, such as masses, radii, and temper- 
ature distributions (Adams et al. 1988; Strom et al. 1989; 
Beckwith et al. 1990). Such spectral modeling, however, 
has revealed an extreme class of T Tauri stars, namely 
flat-spectrum T Tauri stars. The large mid- to far-infrared 
excesses of flat-spectrum T Tauri stars require their disks 
to have a temperature distribution of the form T oc 
(Adams et al. 1988), where R is the distance from the ro- 
tation axis, whereas a T (x R~^/'^ dependence is predicted 
by both a standard accretion disk (Lynden-Bell, Pringle 
1974) and a spatially flat reprocessing disk (Adams, Shu 

1986) . This implies that the disks of flat-spectrum T Tauri 
stars are warmer in the outer regions than predicted by the 
simple disk models; we must therefore consider the mech- 
anisms which heat the outer region of the disk. Disk flar- 
ing allows a disk to receive more emission from the central 
star, and to produce temperature distributions shallower 
than T ex R'^/"^ (Kusaka et al. 1970; Kenyon, Hartmann 

1987) , but does not suffice to reproduce the observed flux. 
The infrared excesses may originate from another cir- 
cumstellar dust component. Calvet et al. (1994) in- 
voked infalling envelopes, which were originally applied 
to the spectral modeling of protostars (Adams, Shu 1986; 
Kenyon et al. 1993), and showed that the infalling en- 



velopes can produce the mid- to far-infrared excesses of 
flat-spectrum T Tauri stars. In fact, observations have 
shown that there exists an extended, disk-like structure of 
radius ~ 1000 AU around a typical flat-spectrum T Tauri 
star, HL Tau (Sargent, Beckwith 1991). Furthermore, 
Hayashi, Ohashi, and Miyama (1993) have revealed an 
infalling motion in the disk-like structure of HL Tau, sug- 
gesting that it is a remnant of an infalling envelope. 

Although the infalling envelope model of Calvet et al. 
(1994) is successful in reproducing the flat spectrum, it 
does not take into account the following two important 
effects. First, as pointed out by Natta (1993), envelopes 
scatter and reprocess the stellar radiation toward the disk, 
and thereby alter the temperature distribution in the disk 
(see also Butner et al. 1994; DAlessio et al. 1997). Second, 
the disk, itself, significantly influences the temperature 
distribution in the envelope. Since these two effects are 
coupled with each other, they cannot be treated sepa- 
rately. Hence, to elucidate the substantial mechanism 
for the flat spectrum, the temperature structure of the 
disk and envelope should be solved simultaneously. In 
the present analysis, we consistently handle the radiative 
energy transport in the disk and envelope by solving the 
angle- and frequency-dependent radiative transfer equa- 
tion in two space dimensions by assuming axisymmetry. 

In this paper, we demonstrate that the mid- to far- 
infrared excesses of flat-spectrum T Tauri stars can orig- 
inate from the disk. Since the envelope which heats the 
disk can be as compact as the disk itself, we henceforth 
call it a halo to distinguish it from an extended infalling 
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envelope. The reflection nebula of HL Tau revealed by the 
Hubble Space Telescope (Stapelfeldt et al. 1995) may be 
regarded as an observational counterpart of the halo. 

2. Model 

The disk-halo model consists of a central star, a cir- 
cumstellar disk, and a halo surrounding both of them. We 
first assume the axial and equatorial symmetry. The cen- 
tral star has luminosity L* and mass M*, and radiates as 
a blackbody with effective temperature T*. We consider 
models with = 0.5 M© and = 4000 K, which are 
typical values for pre-main-sequence stars in the Taurus- 
Auriga molecular cloud. We adopt = 5 Lq for models 
without disk accretion, in which all the energy is radi- 
ated by the central star. Models with disk accretion are 
also examined, for which stellar luminosity is taken to be 
= 2 Lq and the disk, itself, has some intrinsic luminos- 
ity. 

The surface density distribution of the disk is assumed 
to be a power law, 

a(i?)=(7i(i?/AU)-3/2 , (1) 

in the range 1 AU < R< 100 AU. Outside this range, we 
make the surface density go to zero smoothly inward down 
to i?in = 0.1 AU and outward up to i?out = 120 AU. The 
value of (Ti specifies the disk mass. We consider models 
with (Ti = 2 X lO'^ g cm~^, which gives 0.02 Mq for the 
disk mass, and is consistent with the minimum mass solar 
nebula (Hayashi et al. 1985). For the given surface density 
distribution, the two-dimensional density distribution in 
the disk is determined by the hydrostatic equilibrium in 
the vertical direction under the gravity of the central star. 

In addition to the radiation from the central star, vis- 
cous dissipation also heats the disk if it is undergoing ac- 
cretion, as observations of classical T Tauri stars suggest 
(e.g., Bertout et al. 1988). We take into account the effects 
of disk accretion in a simple way as follows. Assuming that 
the mass accretion rate in the disk, Mj, is constant with 
radius, we express the energy generation rate per unit area 
as (Lynden-Bell, Pringle 1974) 
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in the range i?in < R< Rout ■ We further assume that the 
energy generation rate per unit mass does not depend on 
the vertical height, which can be written as 



e{R) = D(R)/a{R) . 



(3) 



The intrinsic disk luminosity is then given by integrating 
equation (2) over the disk area as 
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down to the stellar radius i?* to form a boundary layer, 
where disk material slows down from Keplerian to stellar 
rotation velocity. Radiation from a boundary layer ac- 
counts for the ultraviolet and optical excesses of T Tauri 
stars (Lynden-Bell, Pringle 1974). On the other hand, 
our disk model has an inner radius of i?in = 0.1 AU, which 
is well above the stellar surface. The underlying assump- 
tion for this is that the stellar magnetic field truncates the 
disk at several stellar radii, and accretion from the disk 
onto the central star takes place along the magnetic field 
lines (Bertout et al. 1988). A "bright spot" is then formed 
where the accreting material impacts the stellar surface. 
The luminosity of the radiation from this bright spot can 
be as large as the stellar luminosity, and therefore it will 
also be important as an additional energy source for disk 
heating. In our model, disk material flows at the rate Mc 
from the disk inner radius onto the stellar surface, and the 
amount of radiation energy emitted by the bright spot is 
given by 



GALM^ 1 GM*Mc 

R* 2 i?in 



(5) 



which derives from the assumption that the energy dif- 
ference of accreting material between the total kinetic + 
gravitational energy at i? = i?i„ and the gravitational en- 
ergy at R = R^ is released by radiation from the bright 
spot. The stellar radius, i?*, is defined by =47ri?J(TBT^. 
We do not assume that Mc is equal to Md, i.e., we con- 
sider models with unsteady disk accretion. For simplicity, 
the bright spot is assumed to radiate isotropically as a 
blackbody with an effective temperature of 8000 K. 

The density distribution in the halo is supposed to be 
a power law. 



Ph(r)=pi (r/AU)-P 



(6) 



In a classical picture of accretion disks, the disk extends 



for r > 1 AU, where r is the distance from the central 
star and pi is a reference density. To avoid a singular- 
ity at r = 0, we take ph(r-) = pi [1 + p/2 - p/2{r/ AVf] 
for r < 1 AU. A different assumption for cutting off the 
power-law behavior, for instance, keeping the power law 
down to 0.1 AU, did not change the results significantly. 
Since bipolar outflows may evacuate material along the 
rotation axis, we introduce bipolar holes with an opening 
half- angle of 6'bp. The halo density is reduced slowly and 
smoothly inside the bipolar holes by multiplying a factor 
/bp = 3(6'/0bp)^ -2{9/9bp)^ to the original halo density for 
< ^'bp, where 9 is the polar angle. 

The frequency-dependent absorption and scattering co- 
efficients of dust particles are taken from Miyake and 
Nakagawa (1993), assuming that the size distribution of 
dust particles, n{a), with respect to the radius, a, is ex- 
pressed by n(a) oc a^^'^ for 0.01 /um < a < 1 pm and 
n{a) oc a~^'^ for 1 ^m < a < 10 pm. Miyake and Nakagawa 
(1993) assumed that dust particles are composed of sili- 
cate and water ice. We further assume that the dust va- 
porization temperatures are 1000 K and 100 K for silicate 
and water ice, respectively. For simplicity, we adopt the 
approximation of isotropic scattering. 

Given the density distributions and the radiative and 
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viscous heating sources, we calculate radiative equilibrium 
temperature in the disk and halo which satisfies the con- 
dition of radiative equilibrium, 

4tt Kf^B^(r)di^ = c / Kl^'^E^di^ + e , (7) 
Jo ^0 

where i' is the frequency, nf}'^ is the mass absorption co- 
efficient, Bi, is the Planck function, T is the temperature, 
and El, is the monochromatic radiation energy density. To 
solve equation (7), we adopt the variable Eddington factor 
method in a cylindrical coordinate system {R,Z) (Stone et 
al. 1992), in which we integrate the time-dependent radia- 
tion moment equations forward in time until a stationary 
state is achieved. The radiation moment equations are 
closed by introducing variable Eddington factors, which 
are in turn calculated from the solution of the angle- and 
frequency-dependent radiative transfer equation. For ev- 
ery timestep, the disk is assumed to be in hydrostatic 
equilibrium determined by the temperature distribution 
at the previous timestep. Stationary state solutions ob- 
tained in this way then satisfy the conditions of radiative 
and hydrostatic equilibria simultaneously. We take the 
size of the computational box to be 200 AU x 200 AU 
using a 100 X 100 non-uniform spatial grid, and employ 
203 X 101 angle rays and 101 frequency meshes. The accu- 
racy of our radiative transfer method has been confirmed 
by test calculations described in Stone et al. (1992) and 
in Masunaga, Miyama, and Inutsuka (1998). 

3. Results 

3.1. Standard Case 

The halo density distribution, ph{r) = pi (r/AU)^'^/^; 
is a natural consequence of infalling material with a 
constant mass infall rate, M, if pi is given by pi = 
(47r)-iAf(2GM*)-i/2(AU)-3/2. From the spectral mod- 
eling of protostars in the Taurus-Auriga molecular cloud, 
Kenyon et al. (1993) derived a typical mass infall rate of 
M = 4x 10"^ Mq yr-\ which gives pi = 3x 10"" g cm'^ 
for = 0.5 Mq. Noting this value of pi, we first present 
results for a halo model with p = 3/2, pi = 10"^"' g cm~'^ 
and 6'bp = 60°. In a standard case detailed in this subsec- 
tion, disk accretion is not included, and all of the energy 
arises from the central star with luminosity ~ 5 Lq. 

Figure 1 shows the density and temperature distribu- 
tions in the disk and halo. The temperature gradient is 
shallow along the axis of symmetry where the density is 
very low, while it is steep along the midplane where the 
disk is situated. The resulting temperature distribution is 
far from spherically symmetric, as illustrated by the con- 
tours in the teardrop shape. It can also be seen from figure 
1 that the disk is not isothermal in the vertical direction, 
with a surface temperature higher than the midplane tem- 
perature. The vertical structure at i? = 10 AU is shown 
in figure 2a. The disk surface lies at Z ~ 4 AU, which 
is defined as the position above which the density of the 
disk decreases below that of the halo. As illustrated by 
the solid line in figure 2a, the temperature T begins to 
rise at Z ~ 1 AU, reaches its maximal value at the disk 
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Fig. 1. Density and temperature distributions for the stan- 
dard model having pi = 1.0 X 10~^^ g cm~'^, = 60° and 
p = 3/2. The grey scale and contours represent the den- 
sity and temperature, respectively. The innermost contour 
is r = 100 K, and the contours decrease outward in 10 K 
intervals. 

surface and then decreases. This pronounced temperature 
rise at a disk surface can be explained by the frequency- 
dependence of dust-absorption coefficient. Absorption ef- 
ficiency of dust particles is higher for short-wavelength 
radiation than for long-wavelength radiation. The disk 
surface, which is exposed directly to a high-temperature 
radiation field in the halo, is therefore warmer than the 
disk interior (Chiang, Goldreich 1997). 

In figure 2b the brightness temperature, Tb, at fre- 
quency = 10^^ Hz is shown as a function of Z. We define 
Tb according to I^, = 2hu^ / c^[ex.\){hu / k^Tx,) — 1], where ly 
is the specific intensity, and is calculated assuming that 
the system is viewed pole-on. Also plotted in figure 2b is 
the optical depth, r^, at v ~ 10^^ Hz, which is measured 
from Z = 200 AU toward the midplane. As illustrated by 
the dashed line in figure 2b, Tb is equal to T near the 
midplane where Tj^ 3> 1; as Z increases, Tb also increases 
along with T as long as > 1. When Ti, becomes less than 
unity, Tb begins to depart from T, and keeps an almost 
constant value in the halo. Thus, as expected, we find 
that the observed value of Tb is determined by the value 
of T at ~ 1. Since the halo does not have larger 
than unity, we can define the "photosphere" of the disk as 
the position of = 1 at mid- to far- infrared wavelengths. 
Then, the observed brightness temperature at mid- to far- 
infrared wavelengths traces the photospheric temperature 
of the disk. 

Figure 3 shows radial distributions of the midplane tem- 
perature and the brightness temperature at = lO^'' Hz 
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Fig. 2. Vertical structure at _R = 10 AU for the stan- 
dard model, (a) Temperature T (solid line) and density 
p (dot-dashed line) are plotted as a function of Z. (b) 
Temperature T (solid line), brightness temperature T^, at 
11 = 10^^ Hz (dashed line), and optical depth t„ at v = 10^^ Hz 
(dot-dashed line) measured from Z = 200 AU toward the mid- 
plane are plotted as a function of Z. 
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Fig. 3. Radial distributions of the midplane temperature T 
(solid line) and the brightness temperature Tj, at f = 10^^ Hz 
(dashed line) evaluated at Z = 200 AU. The dotted line rep- 
resents a power law Tb(R) = 350 (fl/AU)-^/^ K. 
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Fig. 4. Spectral energy distributions of the standard model 
at four viewing angles i = 0° (solid line), 30° (dashed line), 
60° (dot-dashed line), and 90° (dotted line). The symbols 
represent observational data taken from Strom et al. (1989) 
for HL Tau (open triangles), DG Tau (solid triangles) and 
RY Tau (open circles). 



evaluated at Z ~ 200 AU. We should be reminded that 
the SED is not determined by the midplane temperature, 
but by the brightness temperature. We find that it can 
be well approximated by Tb(i?) = 350 {R/k\])-^/^ K in 
1 AU < R< 100 AU. This brightness temperature dis- 
tribution is comparable to the temperature distribution 
T{R) = 307 (i?/AU)-o-49 K, which was derived for HL Tau 
through spectral modeling (Beckwith et al. 1990). This 
comparison shows that the disk-halo model can produce 
the large mid- to far-infrared excess of the typical flat- 
spectrum T Tauri star. Moreover, from the result that 
the brightness temperature traces the photospheric tem- 
perature of the disk, it follows that the mid- to far-infrared 
excess originates from the disk. 

In figure 4 the model SEDs are shown at several viewing 
angles. As expected from the results described above, 
the model SEDs are nearly flat in the infrared at viewing 
angles of i = 0° and 30°. At i = 60°, the model SED has a 
weak negative slope in the infrared, with a deep absorption 
feature near A = 10 /im made by silicate and water ice (see 
Miyake, Nakagawa 1993). The water ice feature at 3.1 iim 
is also seen at i = 30° and 60°. In the optical and near- 
infrared, the shape of the model SED varies significantly 
with the viewing angle. The optical depth to the central 
star at V band (A = 0.55 /um) is 3, 18, and 45 at i = 10°, 
30°, and 60°, respectively. Thus, the central star is highly 
obscured when viewed at i = 30° and 60°, and the optical 
flux is due to the scattered hght. 

Observational data taken from Strom et al. (1989) for 
three flat-spectrum T Tauri stars (HL Tau, DG Tau, and 
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Fig. 6. Effects of disk accretion on the SED. The results are 
shown at two viewing angles, i = 0° and 30°. (a) Comparison 
of the steady accretion model (dashed line) with the standard 
model (solid line), (b) Comparison of the unsteady accretion 
model (dashed line) with the standard model (solid line). 




Fig. 5. Closeup view of the density and temperature distri- 
butions to the center for the standard model (upper panel) 
and for the steady accretion model (lower panel). For both 
models, the innermost contour is T = 150 K, and contours 
decrease outward in 10 K intervals. 

RY Tau) are also plotted in figure 4, assuming that the 
distance to these sources is d = f40 pc. Compared with 
the observational data, we find that the model SEDs can 
reproduce the overall spectral shape of these sources, al- 
though the model parameters are not tuned for any spe- 
cific object. 

3.2. Effects of Disk Accretion 

In this subsection we examine the effects of disk accre- 
tion on the temperature distribution and the SED. The 
stellar luminosity is taken to be = 2 Lq , and the disk 
accretion is assumed to contribute to a substantial fraction 
of the total luminosity. We first consider a steady accre- 
tion model, i.e., the accretion rate in the disk is equal to 
that from the disk onto the star. From analyses of the ul- 
traviolet and optical excesses, the accretion rate has been 
estimated to be in the range 5 x 10~*— 5 x 10""^ Mq yr~^ 
for classical T Tauri stars (Bertout et al. 1988). Hence, 



we adopt Md = M^, = 5 x 10^^ A/0 yr^^ for the steady 
accretion model. From equations (4) and (5), one obtains 
accretion luminosities of Ld = 0.18 Lq for the disk and 
Lc = 2.46 Lq for the bright spot, respectively. The total 
luminosity is then 4.65 Lq, which is almost the same as 
that of the standard model discussed in the previous sub- 
section, and thus it is possible to directly compare models 
with and without disk accretion. 

Figure 5 shows a closeup view of the density and tem- 
perature distributions to the center for the standard model 
(upper panel) and for the steady accretion model (lower 
panel). In the standard model, the temperature distri- 
bution is determined entirely by radiative transport of 
energy emitted by the central star, and the temperature 
monotonically increases along the vertical direction in the 
disk. On the other hand, the steady accretion model 
shows a more complicated temperature structure in the 
disk. Inside ^ 10 AU, the temperature is higher at the 
midplane in the steady accretion model than in the stan- 
dard model, indicating that viscous dissipation effectively 
heats the disk in this region. Although the temperature 
decreases along the vertical direction to transport radia- 
tively the energy generated by viscous dissipation, it rises 
again toward the disk surface as in the standard model. 
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Note that the temperature distribution of the steady ac- 
cretion model is quite similar to that of the standard 
model in the halo and the disk outside ~ 10 AU, suggest- 
ing that viscous dissipation is not an important heating 
source in these regions. 

Figure 6 compares the SEDs of models with and with- 
out disk accretion. As discussed in the previous subsec- 
tion, the mid- to far-infrared portion of the SED is de- 
termined mainly by the photospheric temperature of the 
outer disk. It can be seen from figure 6a that the SED 
of the steady accretion model does not differ significantly 
from that of the standard model in the mid- to far-infrared 
region. Therefore, the photospheric temperature of the 
outer disk is determined by the radiation energy from the 
central star and the bright spot, not by viscous dissipation 
in the disk. The radiation energy source of the standard 
model is a 4000 K blackbody with 5 Lq, while that of 
the steady accretion model is a combination of two black- 
bodies with different temperatures. The difference in the 
radiation energy spectra appears in the ultraviolet and 
optical portions of the SED. Compared with the standard 
model, the steady accretion model shows more ultraviolet 
and less optical fluxes. The result shown in figure 6a in- 
dicates that the mid- to far- infrared portion of the SED, 
and hence the photospheric temperature of the disk, is not 
sensitive to the spectrum of the central radiation energy 
source. This is because the inner region of the halo imme- 
diately scatters or absorbs/reradiates radiation from the 
center at longer wavelengths. 

We next consider a model in which accretion takes place 
unsteadily, in the sense that the disk accretion rate is con- 
stant in radius, but it is not equal to the rate at which disk 
material flows from the inner edge onto the central star. 
The possibility that accretion flow in HL Tau is unsteady 
was discussed by Lin et al. (1994). The flattened envelope 
around HL Tau is inferred to accrete dynamically onto the 
disk at a rate of 5 x 10~^ Mq yr~^, while the accretion 
rate from the disk onto the central star should be smaller 
than 7 x 10~^ Mq yr"-'^, because the accretion luminos- 
ity cannot be larger than the total luminosity of HL Tau 
(Lin et al. 1994). In the unsteady accretion model consid- 
ered here, we adopt = 5 x 10~^ Mq yr~^, which gives 
an intrinsic disk luminosity of ~ 1.81 Lq. The accre- 
tion rate from the disk onto the central star is taken to 
be Mc = 2 X 10"'^ Mq yr"\ which is one order magnitude 
smaller than Md, and produces a bright-spot luminosity of 
Lc = 0.99 Lq. The total luminosity then becomes 4.80 Lq, 
which is also comparable to that of the standard model. 

Figure 6b shows SEDs for the standard and unsteady 
accretion models. Comparing these two models, we find 
no significant difference in the mid- to far-infrared por- 
tion of the SED. The unsteady accretion model, however, 
emits more ultraviolet and less optical flux, which may be 
attributed to the difference in the spectra of the central 
radiation energy sources. This feature is similar to that 
seen in the steady accretion model. In the unsteady ac- 
cretion model, however, the central star and bright spot 
emit only about 3 Lq to heat the disk, which would not 
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Fig. 7. Effects of the size of the bipolar lobes on the 
SED. (a) SEDs for a model with pi = 1.0 X 10"" g cm-^, 
= 30°, and p = 3/2. (b) SEDs for a model with 
pi =5.0 X 10-15 g cm-3, ^bp = 30°, and p = 3/2. 

be sufficient to produce a mid- to far-infrared flux as large 
as that of the standard or steady accretion model. This 
deficit of radiation energy is compensated by the accre- 
tion energy generated in the innermost region of the disk. 
In general, most of the accretion energy is released in the 
innermost region of the system where the gravitational 
potential well is deep. In the unsteady accretion model 
we adopt, equation (2) indicates that 80% of the intrinsic 
disk luminosity is generated inside 1.2 AU. The accretion 
energy released in the innermost region of the disk is radi- 
ated at short wavelengths. This emission is subsequently 
scattered or absorbed and reradiated by the halo to heat 
the outer region of the disk, producing the large mid- to 
far-infrared excess, as shown in figure 6b. 

Figure 6b also shows that the unsteady accretion model 
produces more flux at radio wavelengths than the stan- 
dard model. At radio wavelengths, our disk model is op- 
tically thin, and thus the emergent flux is related to the 
temperature at the midplane where disk material is con- 
centrated. Our result indicates that, for a fixed total lumi- 
nosity, disk accretion at a rate as high as 5 x 10~^ Mq yr"^ 
can increase the radio fiux by increasing the midplane 
temperature, although it does not modify the mid- to far- 
infrared portion of the SED. 
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Fig. 8. Dependence of the SED on the reference density 
pi. (a) SEDs for a model with pi = 5.0 X 10~^* g cm~^, 
fbp = 60°, and p = 3/2. (b) SEDs for a model with 
pi = 5.0 X 10-15 g cm-3, = 60°, and p = 3/2. 

3.3. Dependence on the Halo Density Distribution 

As illustrated in the previous subsection, the mid- to 
far-infrared portion of the SED is not sensitive to disk 
accretion. We therefore examine how the flat infrared 
spectrum depends on the assumed halo density distribu- 
tion using models without disk accretion. We first test 
the sensitivity of the SED to the size of the bipolar lobes. 
Figure 7a shows SEDs for a model with the opening half- 
angle of the bipolar holes being 0bp = 30°. The other 
parameters for the halo density distribution are the same 
as those of the standard model, i.e., pi = 10"^'* g cm"'^ 
and p = 3/2. Compared with the standard model hav- 
ing Ohp = 60°, the emergent flux is considerably reduced 
at optical and near-infrared wavelengths when viewed at 
i > 30°, because less optical and near-infrared light can 
escape from the central star through the bipolar holes. 
This means that more energy from the central star is 
used to heat the outer region of the disk. As a result, 
the mid- to far-infrared flux increases slightly, compared 
with the standard model. However, the effect of narrow- 
ing the opening half-angle of the bipolar holes may be 
compensated by reducing pi . Figure 7b shows SEDs for a 
model with a lower halo density: pi = 5 x 10^^^ g cm^"^, 
6'bp = 30°, and p = 3/2. As with the standard model, this 



Fig. 9. Dependence of the SED on the power-law index, p. 
(a) SEDs for a model with p = 1, pi = 1.0 X IQ-^^ g cm"'^, 
and 6»bp = 60°. (b) SEDs for a model with p = 2, 
pi = 1.0 X 10-1* g cm-3, and 9bp = 0°. 

model also produces flat infrared spectra when the view- 
ing angle is 0° < i < 60°. However, some differences are 
found in the SED. The absorption feature near 10 pm is 
shallower due to the reduced halo density, and the slope 
of the SED is more steep from the near-infrared to the 
optical. 

Secondly, we examine how the halo density affects the 
SED. In figures 8a and 8b, model SEDs are shown for 
pi = 5 X 10~^^ and 5 x 10~^^ g cm""^, respectively. For 
both models, we adopt 6'bp = 60° and p = 3/2. The model 
with a higher halo density, pi = 5 x 10^^"* g cm^'^, shows 
relatively narrow SEDs when viewed at i > 0, and the 
peak of the SED occurs in the far-infrared. The emergent 
fiux steeply declines from the near-infrared to the optical, 
because light from the hot inner regions is highly extin- 
guished by the dense halo and only a small amount of scat- 
tered photon can escape the halo. Thus, this model is not 
classified as a "class H" object with fiat infrared spectrum 
unless it is observed precisely pole-on. On the other hand, 
a model with a lower halo density, pi = 5 x 10~^^ g cm"'^, 
still produces flat SEDs at a wide range of viewing angles, 
-30°-60°. 

Finally, we examine the dependence of SEDs on the 
power-law index, p. A molecular cloud core which col- 
lapses to form a star and a disk is expected to be ro- 
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tating. At later stages of the collapse, infalling material 
lands onto the equatorial plane, rather than onto the cen- 
tral star. Thus, the halo density is reduced inside the 
centrifugal radius, and it decreases more slowly with the 
radius than spherically symmetric cases. Wc then consider 
a model with p= 1, pi = 10~^^ g cm~^, and 0bp = 60°. As 
shown in figure 9a, this model also produces flat infrared 
spectra at viewing angles of i ~ 30°— 60°. Compared with 
the standard model, it emits more optical flux, because 
the halo density is lower in the innermost region, so that 
more optical light escapes the halo. 

If the central star drives a stellar wind with constant 
velocity and mass outflow rate, a power law, p oa r^'^, is 
appropriate for the halo density distribution. Figure 9b 
shows SEDs for a model with p = 2, pi = 10~^^ g cm~^, 
and ^bp = 0°. Because the slope of the SEDs is positive 
in the infrared, p = 2 is not favored by the flat spectrum. 

4. Discussion 

4-1. Origin of the Flat Spectrum 

The present results show that the disk-halo model 
can reproduce the overall spectral shape of typical flat- 
spectrum T Tauri stars. In our model, the mid- to far- 
infrared excesses of flat-spectrum T Tauri stars originate 
from the disk, rather than from the halo. The halo needed 
to heat the disk can be as compact as the disk, itself, and 
we suggest that a reflection nebula often associated with 
a T Tauri star is an observational counterpart of the halo. 
The reflection nebula probably corresponds to the inner 
part of a remnant of an infalling envelope. The disk-halo 
model is therefore consistent with the notion that the flat- 
spectrum T Tauri stars are in a transitional stage from 
protostars to T Tauri stars (Hayashi et al. 1993). 

In contrast with the disk-halo model, the infalling enve- 
lope model of Calvet et al. (1994) explains the mid- to far- 
infrared excesses by the emission from infalling envelopes. 
Indeed, observations of molecular line emission have re- 
vealed extended, disk-like envelopes associated with sev- 
eral flat-spectrum T Tauri stars (Hayashi et al. 1993; 
Kitamura et al. 1996a; Momose et al. 1996). However, 
recent high-resolution observations have shown that the 
continuum emission at submillimeter and millimeter wave- 
lengths originates entirely from a compact (^ 100 AU) re- 
gion, suggesting that the extended (~ 1000 AU) envelopes 
traced by molecular line emission do not contribute to the 
submillimeter and millimeter emission (Lay ct al. 1994; 
Kitamura ct al. 1996b; Looncy et al. 2000). Since infrared 
emission originates from warmer regions than submillime- 
ter and millimeter emission, it is quite unlikely that the 
extended envelopes contribute to the infrared emission. 
Hence, the mid- to far-infrared excesses of flat-spectrum 
T Tauri stars should also originate from a compact region. 
The disk-halo model produces emission at any wavelength 
from a compact region, and therefore is more consistent 
with this observational constraint than the infalling enve- 
lope model. 

If the mid- to far-infrared excesses of flat-spectrum 
T Tauri stars are to be explained by the emission arising 



at 100—1000 AU in their infalling envelopes (Calvet et 
al. 1994), they should have been surrounded by sufflcient 
amount of infalling material at a radius of 100—1000 AU 
as protostars are. Growing observational evidence, how- 
ever, indicates that at these radii the amount of infalling 
material around flat-spectrum T Tauri stars is much less 
than that around embedded protostars (Ohashi et al. 
1991, 1996), suggesting that the amount of infalling ma- 
terial inside 100—1000 AU around the central sources de- 
creases during the course of evolution from protostars to 
T Tauri stars. Dispersing motion in the envelopes has ac- 
tually been observed for T Tau (Momose ct al. 1996) and 
DG Tau (Kitamura et al. 1996a). HL Tau is suggested to 
be situated in the bubble wall, which is an expanding shell 
with XZ Tau being the closest known source to the center 
(Welch et al. 2000). In addition, Momose (1998) showed 
that envelopes around flat-spectrum T Tauri stars are less 
centrally condensed than those around protostars. These 
results imply that infalling material around 100—1000 AU 
has been dissipated away, or has accreted onto the central 
star/ disk system. 

In our view, the inner parts of the envelopes around flat- 
spectrum T Tauri stars are still optically thick to stellar 
radiation, although they do not have sufficient mass to 
produce the mid- to far-infrared excesses. For instance, 
the visual extinction toward HL Tau was estimated to be 
Ay > 22 mag (Stapclfeldt et al. 1995), which is comparable 
to that in our calculations. In the disk-halo model, visual 
extinction varies widely depending on the viewing angle 
as well as on the density distribution in the halo; thus, 
our model can explain the flat spectrum for objects with 
various visual extinction values. The detailed structure of 
the halo, however, cannot be predicted by the analysis of 
SEDs, because the halo structure is not unique to give the 
observed flat spectra. Obviously, to solve the degeneracy, 
further observational information is needed on the struc- 
ture within 100 AU. High-resolution observations of the 
scattered light in the near-infrared may provide crucial 
information on the structure of the disk and halo (Close 
et al. 1997). We will report on a comparison between the 
disk-halo model and near-infrared observations in a future 
paper. 

4-2. Effects of Shock Heating 

If the halo is an innermost part of the infalling enve- 
lope, the envelope gas should shock the disk surface. This 
shock heating could be an important energy source which 
our analyses do not take into account. To clarify quantita- 
tively how shock heating modifies the SED in the infrared 
region, detailed radiation hydrodynamics calculations will 
be needed. Although such analysis is beyond the scope of 
this paper, we can discuss that the amount of energy re- 
leased at the disk surface by the shock would not be suffi- 
cient to modify the SED obtained in our standard model. 

Suppose that envelope gas freely falls onto the disk 
along a ballistic orbit, the initial cloud core is spherically 
symmetric and uniformly rotating, the infall is initiated 
by the expansion wave (Shu 1977), and the gas joins a 
Keplerian disk behind the shock. The amount of energy 



Disk-Halo Model for Flat-Spectrum T Tauri Stars 



9 



dissipated at shock per unit mass is GM^,/2R. We denote 
by F{R) the mass flux from the envelope on the disk at ra- 
dius i?, which is given by equation (2.1) in Nakamoto and 
Nakagawa (1994). If we adopt M = 0.9 x IQ-^Mq yr"! 
(Hayashi et al. 1993), i?out = 100 AU (i?out is denoted by 
Rd in Nakamoto. Nakagawa 1994), and = 0.5Mq, the 
energy generation rate due to shock heating at the disk 
surface per unit area, Sg, is given by 

2' 

^^out \ 1 



£s = - - 



R 

1 GM 

2 i?oiit 



■F{R) 



R 



0.45 



R 



100 AU 



■\/l - i?/ i?out 
2 

(l-i?/100 AU) 



-1/2 



-2 -1 

erg cm s 



(8) 

If there were no irradiation from the central star, the disk 
surface would be heated by shock heating to the temper- 
ature 

Teft = (Ss/o-b)^^'' 

- 94 (^^) (1 - i?/100 AU)-i/« K , (9) 

which is much lower than the temperature obtained by 
our standard model including irradiation (see figure 3). 
Therefore, it seems that shock heating at the surface of 
the disk has only a negligible effect on the observed SED 
from the flat-spectrum T Tauri disk. 

5. Conclusions 

1. Wc have shown that disks heated by the scattering 
and reprocessing of the stellar radiation through the 
halo can have flat infrared spectra. 

2. Local viscous heating is not sufficient to produce 
large mid- to far-infrared emission from the disk if 
we consider a reasonable rate of mass accretion in 
disks around classical T Tauri stars. However, the 
accretion energy released in the innermost region 
of the star/disk system is transported by radiation 
through the halo to heat the outer region of the disk, 
resulting again in a flat spectrum. 

3. We examined the sensitivity of the SED to the as- 
sumed halo structure, and found that density distri- 
butions with a power-law index < 3/2 can provide 
the backwarming needed for flat infrared spectra. 
However, we have found that it is difficult to con- 
strain the halo structure only from the SED in the 
infrared. 

4. We have discussed that the halo will be observed as 
a reflection nebula often associated with a T Tauri 
star, indicating that it is the inner part of a remnant 
of an infalling envelope. A more detailed test of the 
model will be made by comparisons with imaging 
observations of near-infrared scattered light. 
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